Suppression of Rabi oscillations for moving atoms 
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The well-known laser-induced Rabi oscillations of a two-level atom are shown to be suppressed 
under certain conditions when the atom is entering a laser-illuminated region. For temporal Rabi 
oscillations the effect has two regimes: classical-like, at intermediate atomic velocities, and quantum 
at low velocities, associated respectively with the formation of incoherent or coherent internal states 
of the atom in the laser region. In the low velocity regime the laser projects the atom onto a pure 
internal state that can be controlled by detuning. Spatial Rabi oscillations are only suppressed in 
this low velocity, quantum regime. 

PACS numbers: 03.65.-w, 32.80.-t, 42.50.-p 



I. INTRODUCTION 

The coupling between the center-of-mass motion of the 
atom and localized laser fields is a central topic in quan- 
tum optics, and leads to many important effects and ap- 
plications in cooling, trapping, deflection, and isotope 
separation experiments. The objective of this paper is 
to point out one further consequence of such a coupling: 
the suppression of Rabi oscillations and the related pos- 
sibility of controlled projection onto internal pure states 
by the motion of the atom from a laser-free region into a 
laser-illuminated region. 

For a general two-level system with an interaction- 
picture Hamiltonian of the form 



H 



(1) 



a state \I> = (*|2)), beginning in the ground state |1) 
(p) at t = 0, will evolve according to 
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where 



(2) 



(3) 



and n is the usual Rabi frequency. 

There is currently much interest in the dynamics of 
two-level systems governed by the Hamiltonian of Eq. 

plus an additional time-dependent driving term, in 
particular to determine the Rabi oscillation suppression 
for critical values of the external field parameters. In 
the so called dynamical localization effect, for example, 
the system remains in one of the two levels . Another 
Rabi-oscillation-suppression effect has been described re- 
cently for the scattering of atoms by a standing light wave 
0. In this case, the partial contributions from different 



diffraction angles lead to different oscillation periods and 
thus to inhomogeneous broadening. 

In this paper we describe a different type of oscillation 
suppression due to the motion of the system rather than 
to an additional time dependent field, or to a diffraction 
effect. Explicitly, we shall consider the motion of a two- 
level atom with transition frequency u) in one dimension, 
with a classical electric field (laser of frequency ujl) illu- 
minating the half axis x > perpendicularly. Without 
damping, and in an interaction picture for the internal 
degrees of freedom, the Hamiltonian can be written in 
the form 



H = f/2m+ - 



ne{x) 

n&{x) -26 



(4) 



where Q is the step function, n plays the role of a laser- 
atom coupling constant, 6 — lul — lu is the detuning, and 
hats denote operators whenever confusion is possible with 
a corresponding c-number. In particular, p is the momen- 
tum operator for the x direction. If one takes damping 
due to photon emissions into account, it can be shown, by 
means of the quantum jump approach "s*], that in three 
dimensions the atomic time development between emis- 
sions may be described by an (effective, non-hermitian) 
"conditional" Hamiltonian He, 



He = p^/2m+-ne{x){\2){l\e''''-y + h.c.} 
- ^{2S + ^J)\2){2\, 



(5) 



where 7 is the Einstein coefficient of level 2, i.e. its 
decay rate or inverse life time, and p is the momen- 
tum operator in three dimensions (3D). The factor e^'^'-^ 
takes into account the spatial dependence of the laser 
coupling. A Hamiltonian of the form of Eq. Q is ob- 
tained from Eq. ((Sj) by neglecting spontaneous emissions, 
i.e. setting 7 = 0, and assuming in addition that the 
atomic wave packet is centered at ?; = and satisfies 
k^Ay 1, at least for some time, so that the exponen- 
tials can be dropped and a one dimensional kinetic term 
suffices. This approximation and its limitations will be 
further commented on Appendix IXI 



2 



The conditional Hamiltonian He is closely related to 
waiting times between photon emissions Indeed, let 
an atomic state |^'(0)) be prepared at time 0. Then it 
can be shown by means of the quantum jump approach 
[3 that 



Po(i) = ||e-''^^(*)/'"'|*(0))|p 



(6) 



gives the probability of no emission until time t and 
n(i) = —Pl^it) is the probability density for the first pho- 
ton. After an emission the atom has to be reset, here to 
its ground state, and then the conditional time develop- 
ment resumes. In this way one can simulate an emission 
sequence and the corresponding Bloch (master) equation 
for the atom. 

As we will show, for moving atoms there are Rabi os- 
cillations both in time as well as in space. In the tem- 
poral case one asks for the probability, P2{t), of find- 
ing the atom in the excited state, without regard to its 
spatial position. For the simple model in Eq. Q with- 
out damping, this can be easily calculated and Figure 
n depicts a suppression for this case. The solid line 
shows the population of the excited state P2 versus t 
for (5 = and for a wave packet which is prepared at 
t = with negligible negative momentum components 
as a minimum-uncertainty-product Gaussian, with the 
atom in the ground state, and far from the laser. At 
sufhciently large times, and for the parameters we have 
chosen, the wave packet is completely within the laser 
region, but there are no Rabi oscillations. Instead, P2 
increases monotonically, and saturates at 0.5; this might 
look surprising at first sight: there is no classical aver- 
aging (the state is quantum and coherent), there is no 
dispersion related broadening, there is no time depen- 
dent term in the Hamiltonian, there is not even decay 
that could smooth or suppress the oscillation. 

Basically, for measurements in time the effect has two 
different regimes and explanations depending on the inci- 
dent energy: at very low kinetic energy, the suppression 
is due to the formation of a pure non-oscillating internal 
state of the atom in the laser region, whereas at interme- 
diate energies it may be understood in terms of a simple 
semiclassical approximation. Yet at higher energies the 
effect disappears, and the expected Rabi oscillations may 
be seen more and more clearly for increasing velocities, 
as illustrated by the dashed and dotted lines of Figure ^ 

Spatial Rabi oscillations can occur in the emission 
probability at different atomic positions. If P2{x,t) de- 
notes the probability density to find the atom in its ex- 
cited state at the position x at time t and if one defines 
I{x) by 



I{x) = / dtP2ix,t), 



(7) 



then "fl(x)dx gives the number of emitted photons (per 
atom) with the atom at dx, and thus I{x) is propor- 
tional to the photon intensity for the atom at position x 
1^ . Neglecting damping, the expression on the rhs can be 




FIG. 1: Probability of excited state P2 versus t — to, 
where to is a mean entrance time, see the definition below. 
fl = 166.5 X lO'' s~\ 7 = 5 = 0. The initial states (at 
t — 0) for the center of mass are minimum-uncertainty 
product Gaussians with A^: = 0.24 /im, (x) = — 1.32^m, 
and (v) — 9.03 m/s (solid line), 49.68 m/s (dotted line) 
and 36.13 m/s (dashed line), to — —{x)/{v). The atomic 
mass is taken to be that of Cs here and in all following figures. 
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FIG. 2: Probability density of the excited state, integrated 
over time, for two different wave packets: (v) = 0.0090 m/s 
(solid line) and (v) = 9.03 m/s (dashed line); n = 166.5 x 10*^ 
s^^; 7 = 5 = 0. The initial state is a minimum-uncertainty- 
product Gaussian in the ground state with = 0.24 ^m, and 
(x) = — 1.32/im. With these parameters Vc = kch/m = 0.28 
m/s and vr = knti/m = 32.3 m/s, see Eqs. 1231 and 1261 . 



calculated with Eq. I^J, and an example is given in Fig. 
13 where the spatial Rabi oscillation is completely absent 
in one of the two cases depicted. In contrast to the tem- 
poral case, spatial oscillations are only suppressed by the 
quantum mechanism, i.e., at very low kinetic energies. 

In the next four sections we shall characterize, neglect- 
ing damping due to spontaneous emissions, several as- 
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pects of the Rabi oscillation suppression due to atomic 
motion: the basic theory (Section^, the effect of the ve- 
locity regimes in the Rabi oscillations for measurements 
in time (Section IIIII) , the effect of detuning and the ve- 
locity in the internal atomic states fSection llV|) . and the 
suppression of spatial Rabi oscillations (Section 0. The 
effect of damping will be taken into account in Section 
ED 



II. BASIC THEORY 

A. Exact results 

Here we shall present the theory required to describe 
the interaction of the moving atom with the perpendic- 
ular laser. We first neglect damping and solve the eigen- 
value equation for the Hamiltonian of Eq. Q subject to 
the condition that the atom impinges on the laser beam 
from the left in the ground state. The eigenvalues are 
E = Ti^k'^ /2m and the eigenfunctions are denoted by 



,(2)j' 



with fc > 0. 



For X < 0, is of the form 



*fc(a;) 



where 



1 



'2tt 



Rie 



x<0, fc > 0, 



(8) 



g2 = fc^ + 2mS/h 



and Im((7) > 0. 
k < \2mS/h\^^^- 



(9) 

For large enough (negative) detuning, 
q becomes purely imaginary and the 
reflected wave for the excited state decays exponentially. 
In the laser region, let \X^) and |A_) be the eigenstates 

of the matrix ^ ^ q 

X±. One easily finds 



corresponding to the eigenvalues 



A± = 
|A±) = 



1 



n 



(10) 
(11) 



Note that |A±) have not been normalized. 

For a; > 0, one can write $fc as a superposition of |A±), 



V2^*fc(x) = C+|A+)e*^'+" + C_|A_)e"'-". (12) 

[The mathematical solutions e^**^*^ are not included be- 
cause they correspond to negative momenta or increasing 
exponentials.] 

From the eigenvalue equation H^k = Ek^k it follows 
that 



(13) 



with Im fc-j- > 0, and from the continuity of $fc(x) at 
X — one obtains, with |1) = (J) and |2) = (°), 

1 + Ri = C7+(1|A+) +C_(1|A_) 
i?2 = C+(2|A+)+C7_(2|A_). 



Similar relations result from the continuity of $5j(x) at 
X = 0, yielding 

C+ = -2k{q + k-)X-/D, 

C- = 2k{q + k+)X+/D, 

Rl = [X+{q + k+)(k-k-)-X-{q + k-){k^k+)]/D, 

R2 ^ k{k^ - k+)n/D, (14) 

where 

D = {k + k^){q + k+)X+ - (fc + fc+)(g + fc_)A_. (15) 
Thus Eq. (|12(l becomes, in components and for a; > 0, 



0i''(^) 



2k 



{{q + fc_)A_e'^+" ~{q + k+)X+e'''-^} 



V2ttD 

kn 



2ttD 



{{q + A:_)e^^-+" ~ {q + fc+)e^'=-"}. (16) 



Time dependent wave packets incident on the laser re- 
gion with positive momentum components and with the 
atom in the ground state may be formed by linear super- 
position 

^{x,t)^ dke-'^''"'/^"'^k{x)4'{k), (17) 
Jo 

where V'(^) is the wave number amplitude at time t = 
corresponding to the freely moving packet. 

The conditional Hamiltonian He for the one- 
dimensional case with damping can be treated in a simi- 
lar way, and this is done in Appendix IbI 



B. Different approximation regimes 

Depending on the incident velocity different approxi- 
mations are applicable. For large velocities a semiclassi- 
cal approximation is sufficient whereas for slow atoms it 
is essential to retain the quantum nature of the transla- 
tional motion. 



1. Fast atoms: semiclassical approximation 

For "high" kinetic energy, fc^ ^ \m{S±n')/h\, we may 
approximate the wave numbers in the wave function by 



q 



mS 



(18) 
(19) 



and therefore, 



/2-K 



,iSt/2 
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2hk 
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in 



ft' 



■ sm 



2hk 



>. (20) 
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FIG. 3: Densities for ground (solid line) and excited state 
(dashed line) ai t — 0.5 /is. The initial state (at f = 0) is 
a minimum-uncertainty product Gaussian for the center of 
mass, with A^; — 0.2436 /xm, (x) = — 1.32/xm, and (v) — 9.03 
m/s; P2{t = 0) = 0; n = 166.5 x lO'^ s'S 7 = (5 = 0, and the 
mass corresponds to a Cs atom. 



[The relation with the Raman-Nath approximation is dis- 
cussed in Appendix El] Up to the normaUzation factor 
this reveals a "spatial Rabi oscillation" that may be as- 
sociated with the temporal one by the change of variable 
t = xm/kh, the time that a classical particle with mo- 
mentum kfi needs to travel from the origin to x. With 
this change the terms in curly brackets correspond ex- 
actly to the internal state amplitudes of the atom at rest, 
see Eq. while the multiplying plane wave denotes an 
undisturbed center of mass motion with fixed momen- 
tum. 

Moreover in this regime reflection is negligible. 



Ri 



R2 



(21) 



A further simplification is to assume that the interfer- 
ences among different momenta are either negligible or 
average out. Then we may quite accurately approximate 
the center of mass motion of a wave packet with high 
kinetic energy by a purely classical free motion, whereas 
the internal state is treated quantum mechanically. More 
specifically, in this approximation the wave packet width 
and spreading are taken into account by using an en- 
semble of semiclassical atoms with their center of mass 
distributed in classical phase space according to the ini- 
tial quantum Wigner function W{qo,po', t = 0). [For free 
motion, the quantum and "classically evolved" Wigner 
functions coincide at all times if initially identical.] To 
each of these atoms we associate an internal state, that 
starts to oscillate once the center of mass crosses x — 
at —qom/po- The expectation values in the semiclassi- 
cal approximation are then computed by averaging over 
the ensemble of semiclassical atoms (translationally clas- 
sical, internally quantum). In particular, the probability 



to find the atom in the excited state at a time t is given 
by 

dqo / dpoW(qo,Po;t^Q) (22) 
-00 Jo 

X (^^^ sm'^[n'{t + qom/po)/2]e{t + qom/pa), 

where we have assumed that i > 0, that the wave packet 
at t = is far from the laser, and that all atoms have 
positive momentum. 



2. Slow atoms 



At the critical wavenumber k = kr 



(23) 



the de Broglie wavelength for the atomic motion 
ArfB(fc) = 27r/fc becomes equal to the "spatial period" 
of the Rabi oscillation in Eq. (|20|l . or "Rabi wavelength" 
A_R(fc) = 27rfc/fc^. For smaller k the quantum aspect of 
the atomic motion cannot be ignored, and k- becomes 
purely imaginary, 



(24) 



whereas k+ remains real. As a consequence, e**^-^ 
evanesces, and the contribution of |A_) beyond k~^ (the 
penetration length for this component) vanishes. 



a 



^e*''■+^ x > k-^ 



/27r 



(25) 



This results in a pure, and non-oscillating internal state, 
as discussed later in Section Hvl 



III. 



TEMPORAL RABI OSCILLATIONS: 
VELOCITY EFFECTS 



In addition to the separation between "classical" and 
"quantum" regimes at fee, see Eq. the classical 

regime may be also subdivided into intermediate and 
high velocities, depending on whether or not the tempo- 
ral Rabi oscillations are suppressed. (We shall see later 
than for spatial oscillations this later subdivision does 
not apply.) 



A. High velocities: (fc) > kn 

Suppose that the center of mass motion of the wave 
packet can be reproduced with an ensemble of classical 
atoms, as discussed in the previous section, and that they 
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FIG. 4; Visibilities, Eq. I|27|l (the quantum and semi- 
classical results are indistinguishable in the scale of the fig- 
ure) for (from left to right) ^ = 0.413, 0.827, 1.654, 2.480, 
and 3.307 xlO^ s'^ . In aU cases the initial Gaussian has 
= 0.2438 /xm, (x) = -1.32 /im, and 7 = 5 = 0. Also 
drawn are the corresponding values of the transition veloci- 
ties vr — knh/m with circles. 



enter into the laser region in a very short time compared 
to the the Rabi period Tr = 27r/i7' (sudden entrance). 
In that case they will oscillate in phase and the Rabi 
oscillations will be seen, as in the dotted line of Figure 
n For Gaussian packets, the transition (average) wave 
number between adiabatic and sudden entrance regimes 
may be identified by equating Tr and a measure of the 
duration of the wave packet passage across the origin, 
tp w 5Axm/ {k)h, where A^; is the wave packet width 
(square root of variance) at the peak's passage across 
X — 0, and the factor 5 is rather arbitrary but is only 
intended to give an estimate, 



2t^ 



(26) 



B. Intermediate velocities: kc < (k) < ka 

In the intuitive language suggested by the classical ap- 
proximation for the center of mass motion, the "mech- 
anism" that explains the suppression of temporal Rabi 
oscillations at intermediate velocities is the averaging of 
the Rabi oscillations of the individual atoms forming the 
ensemble when the entrance time for the whole ensemble 
is long compared to the Rabi period. Since each atom 
enters at a different instant in the laser region, the phase 
of the oscillation will be also different. This is the case 
shown in Fig. O where the peaks of the densities of the 
ground and excited state alternate with a spatial period 
given by h{k)Tfi/m. 

Figure 01 shows the "visibility" Vt obtained from max- 
ima (max) and minima (min) of the probability of the 
excited state, P2{t), corresponding to large times beyond 



the initial transient 



^ P2{max) — P2{min) 
* P2{max) + P2{min) ' 



(27) 



versus the incident average velocity for five different val- 
ues of n and a fixed wave packet width, as well as the 
estimate given by Eq. H26I) : again 6 — 0. The visibilities 
calculated exactly or with the classical approximation are 
indistinguishable. 



C. Low velocity: (fc) < fee 

Below kf. there are no Rabi oscillations either but, ac- 
cording to the discussion of section III B 21 the reason is 
not an effective "averaging" (the scmiclassical approxi- 
mation is not valid), but the fact that the atomic internal 
state formed by the laser does not oscillate at all. 



IV. INTERNAL STATES: EFFECTS OF THE 
VELOCITY AND DETUNING 

The motion of the atom from the laser-free region to 
the laser-illuminated region has very different effects on 
the internal state depending on the incident velocity and 
the detuning. We shall discuss now the effect of the var- 
ious velocity regimes in the degree of mixing of the nor- 
malized, internal, reduced density operator for the atoms 
in the laser region, which is defined by its matrix elements 



J.^dx [^W(x,t)]*^0)(x,t) 
/o°°da;|V'(i)(x,t)|2 + |7/,(2)(x,t)|2 



(28) 



distinguishing between the cases of zero and non-zero de- 
tuning. 



A. Zero detuning 

(i) For (fc) < kc the populations of both ground and 
excited states are equal and the overlap in space of the 
respective wave functions is maximum. After a transient 
time, necessary for Eq. (|25|l to apply, is finally given 
by the pure state 2-^/'^{\l) - |2)), 



P 



1/2 -1/2 
-1/2 1/2 



(29) 



which is nothing but |A-(-)(A-)-|, now normalized. 

(ii) For kc < k < kji excited and ground components 
do not overlap, see Fig. so tends to a diagonal 
matrix. The two populations are still equal (there is no 
global Rabi oscillation because of the semiclassical aver- 
aging effect) and p^ is 1/2 times the unit matrix, which 
amounts to the maximum degree of mixing (incoherence) 
allowed for the reduced internal state. The sharp transi- 
tion from a pure state to a maximally mixed state around 
{k) = kc may be seen in Figure |S1 
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FIG. 5: Degree of mixing versus [v] for 7 = 0, S7 = 3.3 x 10 
s-\ (2;) = -1.34 ^m, and A^; = 0.2436 ^tm. 5 = (solid 
line), — f2/2 (short dashed line) and — f2 (long dashed line). 
Also marked are the corresponding values of Tikii/m (circle, 
square and diamond respectively). 



FIG. 7: Reflection probability to the excited state versus v for 
Q. = 166.5 X 10" s-i; 7 = 0; (5 = n/4 (dashed line), (solid 
line), and —(1/4 (dotted-dashed line), fikc/m is indicated 
with a circle, square and diamond respectively. The peak 
close to 0.2 is at the point where q passes from imaginary to 
real. 




V (m/s) 

FIG. 6: Reflection probability to the ground state versus v for 
Q. = 166.5 X 10® s-^ 7 = 0; 5 = -n/4 (dotted-dashed line), 
(solid line), and f2/4 (dashed line). Tikc/m is indicated with 
a circle, square and diamond respectively. 



There is also quite a dramatic effect of fee on the reflec- 
tion probabilities, which tend to vanish for higher values 
of fc, see Figures Inland [3 for |-Rip and |i?2p respectively. 

{ill) Finally, for k ^ kfj there is no classical averaging 
effect, and the two populations oscillate periodically and 
coherently in time between and 1. The reduced internal 
state in the laser region becomes pure again, as Figure 
13 shows. Note that the transition from intermediate to 
high kinetic energies at kn is rather smooth, in particular 



in comparison with the sharp /cc-transition. 

If the initial atom is prepared in the excited state in- 
stead of the ground state, different scattering eigenfunc- 
tions have to be used. The calculation follows similar 
lines and provides different expressions for reflection am- 
plitudes and coefficients C±. However, the eigenstates 
|A±) remain the same, and it is still true that below 
k ^ kc only k^ is real, so that the same pure state is 
obtained at low kinetic energy. This means that the low- 
energy motion of the atom into a laser illuminated region 
provides a very simple physical mechanism to project any 
initial internal state, pure or mixed, onto the pure state 
2^1/2(11) - |2)) if (5 = 0. The laser illuminated region 
may in fact be moved and the atom be at rest with the 
same effect. 



B. Non-zero detuning 

Detuning may be used as a control knob to obtain 
states with a smaller fraction of excited state. The (un- 
normalized) pure state obtained at low kinetic energy in 
the laser region is given for arbitrary S by 

*=(_4f)- ^^^^ 

As seen in Figure [SJ for the intermediate regime kc < 
(k) < kji, the classical averaging is not so effective for 
non zero detuning, so that the state is less mixed. Note 
also that negative detuning leads to less reflection at very 
low kinetic energies than positive detuning, see Figs. |H| 
and0 and the existence of an additional critical point as- 
sociated with the transition between imaginary and real 
values of q. 
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FIG. 



Spatial visibility, Eq. 



n = 3.307 X 10'' s"\ 7 = 0, and 5 = 0. 



versus v 



for 



V. SUPPRESSION OF SPATIAL RABI 
OSCILLATIONS 

The quantum suppression of Rabi oscillations below kc 
may also be seen in several space dependent quantities. 
They are, however, not sensitive at all to the semiclassi- 
cal mechanism: these quantities oscillate both below and 
above ka defined in Eq H26|l . 

In Figure ISI we show 



\(f>^j^\x)\'^ (max) — \(f)^i^\x)\'^ (min) 

(2) (2) 

Ivk (a;)p(maa;) + ' {x)\'^ (min) 



(31) 



for maxima and minima with respect to x evaluated be- 
yond the penetration length of the A_ -component. One 
may see a rather abrupt jump from to 1 at fc = /cc- 

A similar jump may be observed for I{x) of Eq. Q. 
For wave packets with components below k^, so that k^ ~ 
kc, k- « ike, and for x beyond the transient, small-x 
region, the t integral can be done to obtain I{x) w (v) /v^- 
On the other hand, for k above kc, I{x) becomes 



I{x) 



m 



dk 



■ sm 



^Ix 
2^ 



|(fc|^(i)(0))p. (32) 



These two limiting cases are depicted in Figure |3 



VI. INCLUSION OF DAMPING 

In this section we shall see that the suppression of 
Rabi oscillations can be detected by means of the flu- 
orescence signal of the spontaneously emitted photons. 
Two cases will be considered: a detection measurement 
of the first fluorescence photon, which is described by 
the conditional Hamiltonian, and a measurement where 
all successive photons are also taken into account, which 
requires a master equation (optical Bloch equations). In 



FIG. 9: Fading (not suppression!) of Rabi oscillations for an 
atom at rest: P2 versus t for 7 — 33.3 x lO'' s~^, Q, — 57, 
(5 = 0. The solid line corresponds to the master equation 
and the dashed one to the conditional Hamiltonian ("first 
photon" ) . 



both cases the fluorescence signal is proportional to the 
population of the excited state P2 which can be calcu- 
lated by means of He or by the Bloch equations, respec- 
tively. 

The suppression of Rabi oscillations discussed here 
should be distinguished from their fading away due to 
damping and approach to a stationary internal state. 
The latter is quite well-known for atoms at rest (cf. 
Figure El . It is possible to operationally distinguish fad- 
ing and suppression by a combination of U, and 7 that, for 
the atom at rest, allows to observe several oscillations be- 
fore reaching the asymptotic stationary population. This 
corresponds to strong driving conditions, f7 3> 7. 



A. Rabi oscillations in the first-photon distribution 

For the one-dimensional case the conditional Hamilto- 
nian He can be written as 



Hr^ 



P 

2m 



ne{x)e-'''^y 

ne{x)e^''^y -2(5-^7 



(33) 



The generalization of Eqs. pill5|) for 7 7^ is provided 
in Appendix IbI 

For high velocities the entrance of the wave packet 
in the laser region is sudden, so that P2 versus time 
takes the same form as for the atom at rest from the 
entrance time on. At lower velocities the oscillation pat- 
tern disappears; kji still marks the transition, which is 
now smoother than when 7 = 0. Figure [1111 shows P2 for 
a velocity below kn. The Rabi oscillation suppression is 
evident, compare with the dashed line of Figure |5| 

The transition at kc is illustrated in Fig. ^] which 
shows the degree of mixing versus the average velocity 
Note the smoothing effect of a non-zero 7. 
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FIG. 10: P2 versus f for 7 = 33.3 x 10*' s-\ = 57 and 
(5 = 0. The parameters of the initial Gaussian are (v) = 9, 03 
m/s, Ax = 0.2436 fim, and (a;) — —1.322 fim. 
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FIG. 11: Degree of mixing of the internal state in the laser 
region versus (i;) for A.x= 0.2436 fim, 7 = 33.3 x 10^ s~^, Q = 
IO7 (solid line), and Q — 5^ (dashed line). The corresponding 
values of Tikc/m are also shown with a circle and a square 
respectively. 



An observation of the suppression of Rabi oscillations 
using the first photon could be achieved for a "Lambda" 
configuration of three atomic levels such that the laser 
couples two of them, while a third one acts as a sink for 
the excited state j^l (only one photon may be emitted per 
atom) . The ideal experiment would require a preparation 
of the atoms, sent one by one, according to a single wave 
function, and a detector capable of responding to a single 
photon. The experiment would be repeated many times 
until a profile similar to the one in Fig. 1101 is obtained 
for the density of photon detection times. 



B. Temporal Rabi oscillations in the photon 
intensity 

For a simple two-level atomic configuration, without 
a sink level for the excited state, the atom may emit 
many photons. In this case recoil effects may become 
important, especially at low incident energies, and lead to 
spatial broadening. The master equation that describes 
the density operator p^^ of the atomic system in three 
spatial dimensions is given, following IQ], by 

= -'-{H,p^°-p^^Hl} + n{p^% (34) 

where He is the conditional 3D Hamiltonian, Eq. (O , and 
the "reset" or "jump" term TZ{p^^) takes into account the 
atomic recoil along the directions weighted by the dipole 
emission distribution. In momentum representation, 

(p|7^(p3^)|p') ^ 7|1>(1| (35) 
X J dKP{K){p + MlkIpI^Ip' + ^kLK.), 

where the integral is over all possible photon directions 
represented by the unit vectors k = k/|k|, and 

3/ K-d\ 3 9 

9 being the angle between k and the dipole moment d, 
which is taken along the z direction. 

Even though the fluorescence photons and the associ- 
ated atom recoil may go in any direction, one can take 
the trace over momentum components Py and Pz in the 
reset term, and reduce the integral to a one dimensional 
one, see e.g. Illj, to obtain for the initial atomic- motion 
direction x the expression 

TZAp) ^ Tr,,,7^(p^^) = |l)(l|7^ (37) 
X J du{l + u'^){px + hkLu\p22\p'x + fikLu), 

where the reduced atomic density operator for the a;- 
direction is 

p,, =Tr,,,(z|p3^|j), i = l,2, J = 1,2. (38) 

Note that it contains both internal and center of mass 
information of the atom, so it is different from p^ in 
Eq. H28|l . Even though, taking the trace over y and z 
in the other terms of Eq. (|34|l does not lead to a closed 
equation for the reduced one-dimensional density matrix 
of Eq. H38() . a one-dimensional approximation is valid, 

p = -^{H,p-pHl}+nx{p), (39) 

with He given by Eq. (|33|) . provided that the effect of 
the two broadening mechanisms in the y direction, the 
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FIG. 12: P2 versus time for 7 = 33.3 x lO'^ s"\ t7 = 57, 
Ax= 0.12 yLxn, (x) — — 1.322/im, {v) = 3.61 m/s; average 
for 10000 trajectories (dashed line), compared to P2 for a 
very fast packet (solid line). The later curve corresponds to 
the solid line of Figure |5] displaced so that the starting time 
coincides with the entrance of the center of the wave packet 
in the laser region. 

standard quantum mechanical wave packet spreading and 
the atomic recoil, remain small. The time scales of both 
effects are estimated in Appendix 1X1 

The fluorescence signal is proportional to P2, which is 
now the total population of excited atoms, regardless of 
whether or not they have emitted. The one dimensional 
master equation may be solved in principle directly, or, 
as done here, using the quantum jump technique, i.e., 
averaging over many "trajectories" For each of them 
the photon detections occur at random instants. Figure 
m shows P2 averaged over 10000 trajectories, compared 
with the solid line of Figure 1^1 properly shifted, which 
gives the time evolution of P2 for a sudden entrance of a 
very fast wave packet. The suppression of Rabi oscilla- 
tions is evident. 

The ideal fluorescence measurement would be per- 
formed by sending one atom at a time, as described be- 
fore. In a less demanding experiment, one could prepare 
a cloud of many atoms of sufficiently low density, and 
send it towards the laser, measuring the fluorescence sig- 
nal versus time. The cloud ensemble would lead to sup- 
pression of Rabi oscillations at sufficiently low velocities. 
In this case the incident translational state would not be 
pure but a true mixture. 

The transition at fc^ would be most easily noticed by 
the reflection dominant at lower kinetic energies. For the 
Cs transition we are considering in the numerical exam- 
ples Vc = hkc/m — 0.28 m/s, which is well above the 
recoil velocity limit hkL/m = 0.35 cm/s. The possibility 
to control the pure state formed in the laser region will 
depend on the ability to focus the atomic beam along 
the y- direction in the scale of the laser wavelength, due 
to the y-dependence of the pure state obtained, see Eq. 
(lB5l . 



C. Spatial Rabi oscillations 

With damping included, P2{x,t) becomes {x\p22{t)\x) 
so that I{x) in Eq. iQ becomes 

/oo 
dt {x\p22{t)\x) , (40) 
-00 

and "fl{x) is the mean number of photons per atom per 
unit length. 

In a similar way, one may consider the probability den- 
sity, denoted by 7/0(0:), for the atomic position when the 
first photon is emitted. Similar to Eq. Q, Iq{x) is given 

by 

/CO 
dt\i:('\x,t)\^ (41) 
-00 

where the time development is given by the conditional 
Hamiltonian in Eq. and ip'^^^{x,t) denotes the 

excited-state component. Integrating Eq. (|41|l over x 
gives the fraction of atoms which emit a photon, so that 
/ dxIo{x) tends to one for kc- 

Above kc, and using the language of the classical ap- 
proximation, the time when the atom arrives at the laser 
and starts the Rabi cycle is unimportant for an observa- 
tion of the spatial dependence of the photon intensity; the 
only relevant information is the position where the emis- 
sion takes place. In other words, the suppression will not 
be visible at the intermediate velocities below kji; for an 
example see Figure IT^ where both 7/0(2:) and jl{x) are 
depicted for the same initial conditions, with no suppres- 
sion of the oscillations. In the opposite "quantum" case, 
k < kc, the suppression effect is visible in Io{x); for an 
example see the dashed line of Figure E| Note the clear 
distinction between the small- a; region, where there is still 
a contribution from |A_), and the larger-x region, which 
depends only on |A+). A similar result holds for I{x) but, 
because of the continuation of pumping-emission cycles, 
there is no exponential decay; see Figure 1141 for an ex- 
ample. We may estimate I{x) after the transient peak 
by approximating {x\p22{t)\x) in Eq. (|40() by its value 
at 7 = 0, as in Section This gives, with k ^ kc, 
I{x) « {v)/vc- The exponential decay of /o(a;) may also 
be approximated in the strong driving limit and k ^ kc 
by retaining the dominant exponentially decaying terms 
in Eq. 123), loix) « (?;)e-^('"/"^)'^'/Vi;2. 

VII. DISCUSSION 

The possibility to reach ultra cold temperatures for 
atomic motion leads to a number of new interesting quan- 
tum phenomena that would remain hidden otherwise. As 
is well known, the internal level populations of fast atoms 
exhibit Rabi oscillations when traversing a sufficiently 
strong traveling laser field. In this paper, however, we 
have shown that for cold atoms these Rabi oscillations 
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FIG. 13: 7/(2;) (computed with 160 trajectories, solid line) 
and 7/0 (dashed line) versus x for 7 = 33.3 X 10'^ s'^ , fl = 57, 
5 = 0, 0.12 fim, {x) = -1.322/im, and {«) = 3.613 m/s 

{vr — kRh/m — 16.15 m/s, and Uc = kcTi/m — 0.28 m/s). 




0,05 

x(|a.m) 



FIG. 14: 7/(2;) (for 40 trajectories, solid line) and 'yloix) 
(dashed line) versus x, same conditions as previous figure ex- 
cept {v) = 0.18 m/s. 



kinetic energy the suppression of Rabi oscillations is as- 
sociated with the projection onto specific internal pure 
states, which can be controlled via detuning. Purifica- 
tion procedures are essential in quantum data processing 
[Tk] . Our results are also relevant for the analysis of time 
scales defining the particle's traversal through a region 
of space ^. They suggest that the excited state popula- 
tion cannot in general be taken as a measure of the time 
elapsed in the laser region (based solely on a comparison 
with the Rabi oscillation of the atom at rest) unless the 
particle's kinetic energy is sufficiently high. 

Even though we have centered the discussion on an 
atom-laser system, the formalism is equally valid for 
other two-level systems, and in particular for a localized 
magnetic field acting on a spin-1/2 particle, i. e., to the 
so called "Larmor-clock" , which is usually applied with 
an additional potential barrier and in the weak coupling 
limit 

A simplification of our treatment has been the semi- 
infinite laser-illuminated region. This is not a major ob- 
stacle to observe the effect since, due to the low velocities 
implied, actual finite laser profiles would be enough for 
the occurrence of many Rabi oscillations of a classically 
moving atom. For example, a Rabi period for our Cs 
transition takes 0.036 /xs when f2 = 57, whereas an atom 
with speed ten times Vc would only move 0.1 /im in that 
period. 
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APPENDIX A: LIMITS OF THE ID MODEL 



may be suppressed and in fact can be totally absent. In 
the context of moving two-level atoms, two types of Rabi 
oscillations can occur, temporal and spatial ones. In the 
temporal case one observes the population of the upper 
level in time and the associated photon emission, regard- 
less of the atomic position. In the spatial case one ob- 
serves where the emissions occur in space, without regard 
to time. 

We have distinguished two regimes for the suppression 
of temporal Rabi oscillations at low and intermediate 
kinetic energies and have characterized the transitions 
quantitatively. Spatial oscillations are only suppressed 
at low kinetic energies. 

Quite generally the effects of slow atomic motion or 
laser displacement on a two-level qubit are of interest in 
the design on quantum computers p^ . Applications of 
the present effect may be based on the fact that at low 



A simple random, classical model will enable us to esti- 
mate a maximum time of validity of the one dimensional 
approximation used in the text. At variance with other 
random walk approaches dealing with stimulated emis- 
sion and absorption processes in standing waves |14| , the 
present model focuses on the momentum kicks due to 
spontaneous photon emissions. 

Suppose that the atom suffers random kicks at j — 
1, 2, 3, n instants separated by the interval Ai and that 
at t = the atom is at rest, y — v — where v is the 
velocity in the y-direction. If Auj and Ay^ are velocity 
and position increments in the y-direction in step j and 
Vy and yj the corresponding velocity and position, we 
have 

y„ = ^Ay, = At^i;^ (Al) 
3 3 



11 



and 

3 



vi = (Ai;i) + (Aui + Au2) + ... + (Awi + 
nAvi + (n - l)Aw2 + {n - 2)Av3 + 



+ AVn) 

+ Aw„, 



so that 



yl = (A02[n2(Az;i)' + (ri-l)'(A«2)' + ... + (At;„)2] 
+ crossed terms involving AviAvi, i ^ I. 

When taking the average, the contribution from the 
crossed terms vanishes, 

(yl) = (A^)2[n2 + (n-l)2 + ... + l]((A^;)2) 



n(n + l)(2n+ 1) 
6 



(AtrHAv)^ 



where ((Aw)^) is the variance of the velocity increments 
in the y direction at any of the j instants. For large n, 

,3 



{y')^^iAtn{Avr). 



(A2) 



For the dipole radiation distribution of Eq. . and 
considering that the modulus of each velocity increment 
is due to the recoil velocity Av = kLfi/m, we find 
{{Avyf) = lAv^jh. A critical n may be obtained from 
(|A2(I by imposing, say, A:i(y^)^/^ — 1/10. For the transi- 
tion of Cs atoms we are considering, 7 = 33.3 x 10® s~^, 
fc^^ = 852/27r nm, and for strong driving At w 2/7; this 
gives n k, 30, or approximately 2 [is of time from the first 
photon emission. 

The other phenomenon to take into account is quan- 
tum mechanical dispersion of the wave packet. This may 
be estimated by assuming 

A;i(y2)i/2 ^ 1/10 in the for- 
mula 



{Ayf = {Ay)l 



2.2 



1 + 



fiH 



4m2(Ay)4 
For kL{Ay)f) = 1/20, this gives 3 fis. 



APPENDIX B: THEORY WITH DECAY 



(A3) 



We provide generalizations of Eqs. H11I15|1 for 7 7^ 
and an arbitrary value of y corresponding to the Hamil- 
tonian of Eq. ifSSj) : 



A+ = 



1/2- 



|A±) 



r 
4 



+ is-f + n'^ 



1 



(Bl) 



(B2) 



^2 ^ ^2 _ 2mA± ^2 _ ,„2 



2m , 



^ , 9 =r + ^('5 + J7/2), (B3) 



where fc± and q have positive imaginary parts. In terms 
of these variables, the expressions for C±, and D are 
the same as in Eqs. H14|l and H15|l . whereas 



R2^k{k+-k^)ne'''''y/D, 



(B4) 



In this one-dimensional approximation the parametric 
dependence on y does not affect the reflection probabil- 
ities, A±, q, k±, or C'±, but it gives a phase factor to 
i?2, and to (2|A±). A consequence of the later is the for- 
mation of different pure states in the laser region at low 
kinetic energy. In particular, for 7 = 0, 



IA+) = 



1 



(B5) 



APPENDIX C: RELATION WITH THE 
RAMAN-NATH APPROXIMATION 

It is interesting to compare the approximation of the 
text (x as a variable, y as a parameter), with the Raman- 
Nath (RN) approximation where y is the variable, and 
X a parameter, expressed in terms of time through an 
assumed (classical) linear relation 



t — mx/kh, 



(CI) 



where k is also a parameter. Consider the atom evolving 
with the Hamiltonian 



^«^-2Ue^^--^ -25 



(C2) 



Note the absence of a kinetic term, a basic feature of 
the RN approximation. This Hamiltonian is easily di- 
agonalized. The eigenvalues are given by Eq. (|10|1 and 
the eigenvectors take the form of Eq. l|B2p so that the 
wave vector for a state initially in the ground state with 
amplitude ipoiu) is given by 



^flAr(y) = My) 



n' -6 



^RNiy) = My) 



2n' 



2VL' 



2^1' 



or "ifn.N = V'o^flN- Making the substitution of Eq. (|C1|) 
in one finds the same result as $k of Eq. 112() when 
using the "high kinetic energy" approximations of Eqs. 
I|18I19|I . except for the plane wave factor e*'^'^/(27r)^/2^ 



(27r)i/2 



RN- 



(C3) 



Physically, the plane wave accounts for the undisturbed 
translational motion along the incident direction whereas 
the internal one is described by the RN wave vector ^b,n- 
In the RN approximation the probability to find the 
atom in the excited state is given by fi^ sin^(tri'/2)/il'2. 
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In spite of the neglect of the spatial displacement along 
y, there is a possibility of momentum exchange due to 
the absorption of a laser photon . Taking the Fourier 
transform we find the following normalized momentum 
distributions for ground and excited state 

Vl{Py) = |(p,|V'o)P, (C4) 

V2[Py) = \{py - kLh\^^)\\ (C5) 



Note that the RN approximation cannot describe the 
Rabi-oscillation-suppression effect at low kinetic energies 
or the transition region around kc since the longitudinal 
motion is treated classically. The suppression effect at in- 
termediate energies could be described if an ensemble of 
"RN atoms" mimicking a wave packet distributed along 
the longitudinal direction is considered. 
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